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1.  INTRODUCTION 

The  Tactical  Environment  Multiple  Systems  Evaluation  Prograun 
(TEMSEP)  is  jointly  funded  by  the  Army  and  the  Defense  Nuclear  Agency 
and  is  directed  toward  vulnerability  assessments  auid  hardening  of 
tactical  military  equipment  for  tactical  nuclear  threats.  TEMSEP 
complements  the  Multiple  Systems  Evaluation  Program  (MSEP) , vdiich  is 
also  directed  toward  vulnerability  assessments  and  hardening  of 
equipment  for  nuclear  threats.  Unlike  TEMSEP,  however,  MSEP  is  directed 
towards  stategic,  rather  than  tactical,  nuclear  threats. 

In  the  absence  of  threat-re latable  simulation  of  source-region  EMP 
environments,  meaningful  assessments  of  source-region  EMP 
vulnerabilities  of  tactical  equipment  are  feasible  only  through 
theoretical  and  controlled  experimental  evaluations  of  dominant 
source-region  coupling  mechanisms.  The  experimental  evaluations  eu:e 
being  performed  at  flash  x-ray  facilities,  such  as  AURORA  and  HIFX  (the 
High- Intensity  Flash  X-Ray) , both  at  Harry  Diamond  Laboratories  (HDD . 
These  facilities  are  being  used  as  partial  simulation  tools  ^uld  as  tools 
for  validating  theory. 

During  the  experiments  at  AURORA  and  HIFX,  parameters  such  as 
coupled  currents,  magnetic  fields,  and  Compton  currents  are  measured 
with  sensors  such  as  current  probes,  Moebius  loops,  and  Rogowski  coils. 
Typically,  the  sensors  are  part  of  a sensor-cable-balun  system  such  as 
the  one  shown  in  the  block  diagreun  of  figure  1. 

Because  it  is  less  noisy  in  a radiation  environment  than  most  other 
cables,  twineix  cable  is  normally  used  to  tr^ulsmit  the  sensor  output  from 
the  test  cell  into  the  data  room.  Unfortunately,  twinax  cable  is 
relatively  lossy,  and  the  signal  loss  is  frequency  dependent.  When 
twineix  cable  is  used,  a balun  is  needed  to  match  the  78-n  impedance  of 
this  cable  to  the  50-ft  impedance  of  the  cable  used  in  the  data  room. 


Figure  1.  Block  diagram  of  measuring  system  used  at  AURORA. 


I 
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Each  sensor  exhibits  a particular  response  to  a given  signal , and 
this  response  is  degraded  by  the  cables  and  balun  through  which  the 
signal  passes.  Thus,  the  voltages  measured  by  the  oscilloscopes  in  the 
data  room  represent  the  response  of  the  sensor- cable-balun  system  to  a 
particular  phenomenon  occurring  in  the  test  cell,  and  these  voltages 
must  be  somehow  "transformed"  into  the  pareuneter  being  measured. 

The  response  of  the  sensor-cable-balun  system  to  input  signals  can 
be  measured  with  a network  analyzer.  This  response,  represented  by  the 
transfer  function  measured  by  the  network  analyzer,  provides  the  ampli- 
tude of  the  output  voltage  cind  the  phase  difference  between  this  voltage 
and  the  input  signal  for  an  input  signal  of  unit  amplitude  over  a wide 
range  of  frequencies.  If  a current  probe  is  used  to  measure  it,  the 
transfer  function  is  actually  the  admittance  of  the  measuring  system, 
cind  the  input  signal  can  be  found  from  the  relationship 


Kw)  = Y(w)V{u)) 


where  I is  the  input  current,  Y is  the  admittance  (transfer  function)  of 
the  measuring  system,  V is  the  measured  voltage,  and  ui  is  the  arngular 
frequency. 

The  computer  program  TRANS2  uses  this  transfer  function  information 
to  transform  oscilloscope  voltage  measurements  into  the  parameters  being 
measured,  even  if  they  are  not  simply  currents. 


2.  DATA-PROCESSING  PROCEDURE 

The  data-processing  procedure  which  is  used  to  accomplish  this 
transformation  is  this:  The  oscilloscope  pictures  are  digitized  to 
provide  information  which  can  be  manipulated  by  a digital  computer.  The 
digitized  information  is  stored  on  a computer  disk  file  and  then 
plotted.  The  plot  is  compared  with  the  oscilloscope  picture,  and  any 
needed  corrections  are  made  to  the  digitized  data.  The  digitized 
information  is  then  read  off  the  disk  file  by  TRANS2,  which 

(a)  multiplies  the  data  by  the  appropriate  sensor  response  and 
conversion  factors, 

(b)  numerically  transforms  the  result  into  the  frequency  domain, 

(c)  multiplies  the  frequency-domain  representation  by  the  appro- 
priate transfer  function, 

(d)  transforms  this  information  back  into  the  time  domain. 
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(e)  integrates  the  result,  if  necessary,  and 

(f)  plots  the  resultant  information. 


The  plotted  information  represents  the  time  histories  of  the  actual 
fields,  currents,  etc.,  occurring  at  the  sensor  locations  during  the 
AURORA  or  HIFX  tests,  and  not  merely  the  measured  voltages.  These  plots 
can  then  be  compaired  with  theoretical  predictions  of  test  results. 


3.  PROGRAM  DFILE 

A listing  of  the  program  DFILE,  which  is  used  to  put  the  digitized 
data  onto  a computer  disk  file,  is  presented  in  appendix  A.  Figure  2 
shows  a sample  deck  needed  to  run  the  program. 


//JFDIADFI  JOB  (HK3020) , MFWDIETZ' ,CLASS=A,MSGLEVEL=(2,I) 
Z/STEPl  EXEC  TRANS, NAME='&6DATA' 

//SYSIN  DO  * 
eOUTPUT  SEND 

data  set 


data  set  llt\ 

(7/8/9)  CARD 
/■■■ 

//STEPI  EXEC  ANAF0RT,0UT=X 

//SYSIN  DD  * 

source  deck 


//G0.FTI2F00I  DD  DSN=* . STEP  I . STP . FTl 2F00 1 , D I SP= (OLD  DELETE) 
//  DCB=(RECFM=VS) 

//GO.FTlAFOOl  DD  D I SP= (NEW, CATLG) . V0L=SER=USER02 . 

//  UNIT=SYSDA,SPACE=(TRK.(2,2)) , 

//  DSN=HK3020.CR362022 

// 

Figure  2.  Sample  deck  for  DFILE. 
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3.1 


Description  of  DFILE 


The  procedure  TRANS,*  which  transforms  the  digitizer  output 
into  usable  data,  is  executed  and  the  data  are  written  to  logical  unit 
No.  12.  DFILE  calls  the  subroutine  READ IT,  which  reads  the  data  from 
logical  unit  No.  12.  The  subroutine  CSTOUT,  contained  in  the 
ANAPAC^  library  and  attached  by  the  ANAFORT*  procedure,  is  then  called. 
CSTOUT  checks  the  time  ordering  of  the  independent  array  and  casts  out 
those  points  not  in  cin  ascending  time  order. 

DFILE  converts  the  data  to  double  precision  to  increase  the 
accuracy  of  the  data-processing  effort  and  calls  WRITED,  which  writes 
the  double-precision  data  to  logical  unit  No.  14  so  that  it  can  be 
cataloged  during  the  GO  step.  Finally,  DFILE  prints  the 
double-precision  data. 

3 . 2 DFILE  Sample  Deck 

A few  comments  concerning  the  Scunple  deck  shown  in  figure  2 are 
now  appropriate.  The  namelist  entitled  "OUTPUT"  contains  two  logical 
variables,  LIST  and  PLOT,  which  are  input  to  TRANS.*  A value  of  TRUE 
for  either  LIST  or  PLOT  turns  on  either  the  printing  or  the  plotting  of 
the  data  sets,  respectively.  This  printing  or  plotting  is  performed  by 
TRANS  and  occurs  before  CSTOUT  corrects  the  time  ordering.  Second,  the 
last  card  in  each  data  set  must  be  a (7/8/9)  card,  and  there  must  be  an 
extra  (7/8/9)  card  after  the  last  data  set.  Finally,  the  name  of  the 
data  file  cataloged  by  DFILE  is  given  by  the  DSN  parameter  on  the  last 
executable  card. 


4.  PROGRAM  TRANS2 

After  the  digitized  data  are  placed  on  file  and  verified  to  be 
correct,  TRANS2  must  be  run  to  transform  the  data  into  the  measured 
fields,  currents,  etc.  TRANS2  is  currently  maintained  as  an  UPDTE^  file 
so  that  temporary  changes  in  the  program  can  be  readily  made.  A seunple 
deck  needed  to  create  an  UPDTE^  file  containing  TRANS2  is  shown  in 
figure  3.  A listing  of  TRANS2  is  presented  in  appendix  B.  Definitions 
of  the  input  parameters  for  TRANS2  are  presented  in  section  4.1, 
followed  by  a description  of  how  the  program  works  and,  finally,  by  a 
discussion  of  the  sample  deck  needed  to  run  TRANS2  (sect.  4.3). 

^Thomas  V,  Noon,  User's  Manual  for  the  Modular  Analysis-Package 
Libraries  ANAPAC  and  TRANL,  Harry  Diamond  Laboratories  TR-1782  (November 
1976) . 

^IBM  OS/VS  Utilities  Manual,  GC35-0005-4  (September  1976). 

*From  recommendations  by  Egon  Marx  and  Thomas  V.  Noon,  Harry  Diamond 
Laboratories  (1976). 


//JFD1CTR2  JOB  (HK3020) , ' JFWDI ETZ ' ,CLASS=C 
//  EXEC  PGM=I EBUPDTE,PARM=NEW 

//SYSPRINT  DD  SYSOUT=A 

//SYSUT2  DD  UNIT=SYSDA,DISP=(NEW,CATLG) ,V0L=SER=USER02, 
//  DSNAME=HK3020.TRANS2.SPACE=(TRK, (25,2)) , 

//  DCB=(RECFM=FB,LRECL=80,BLKSIZE='4080) 

//SYS  IN  DD  * 

./  ADD  LIST=ALL 

./  NUMBER  NEW1=1000, INCR=1000 


source  deck 
./  ENDUP 
/* 

1 1 


Figure  3.  Sample  deck  to  create  UPDTE  file  for  TRANS2. 


4. 1 Input  Parameters  for  TRANS2 
Card  No.  1:  FORMAT  (2110) 


NDSETS  Number  of  data  sets  to  be  processed. 

ISAME  If  equal  to  one,  the  input  parameters  for  the  first 
data  set  are  used  for  all  the  data  sets;  otherwise, 
input  parameters  must  be  provided  for  each  data  set. 

Card  No.  2;  FORMAT  (7110) 

NOFF  Number  of  data  points  removed  from  the  front  of  the 

pulse . 

NOFFR  Number  of  data  points  removed  from  the  rear  of  the 
pulse. 

INVERT  If  equal  to  one,  the  pulse  is  inverted;  otherwise, 
not. 

IZERO  If  equal  to  one,  the  amplitudes  of  the  first  and  last 
data  points  are  set  equal  to  zero;  otherwise,  not. 

IPLOT  If  equal  to  one,  the  input  digitized  data  are  plotted 
after  any  necessary  time  shift  is  made;  otherwise, 
not. 

If  equal  to  one,  the  processed  data  are  written  to 
logical  unit  No.  9 so  that  they  can  be  cataloged  as  a 
disk  file;  otherwise,  not. 
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IFILE 


T 


ISENSR  Denotes  the  type  of  sensor  used  to  produce  the  data: 

1 = Current  probe  used  to  measure  current  on 

truncated  cylinder  banjo. ^ 

2 = Current  probe  used  to  measure  current  on  single 

wire. 

3 = Current  probe  (Rogowski  coil)  used  to  measure 

Compton  current. 

4 = Moebius  loop. 

q 

5 = Current  probe  on  pie  pan. 

P = One-turn  Rogowski  coil. 

Card  No.  3:  FORMAT  (2110) 

IPLTC  If  equal  to  one,  coupled-current  output  from 
SAPSC**  is  plotted  on  the  same  graph  as  the  processed 
data;  otherwise,  not.  Set  equal  to  one  only  if 
ISENSR=1. 

NCUR  If  equal  to  one,  SAPSC  output  is  plotted  for  the 

first  measurement  position;  otherwise,  for  the  second 


measurement  position.  Set  equal  to  one 

ISENSR=1. 

only  if 

Card  No. 

4:  FORMAT  (110,  E10.3,  110) 

NSTAR 

Number  of  frequency  and  time  points  used  by 
Fourier  transform  (FFT)  (must  be  a multiple 

the  fast 
of  two) . 

FMAX 

Maximum  frequency  used  by  the  FFT. 

KPLOT 

If  equal  to  one,  frequency-domain  results  are 

plotted;  otherwise,  not. 

Card  No. 

5:  FORMAT  (3110) 

NTRANS 

Number  of  transfer  function  amplitude  points 
less) . 

(25  or 

^John  F.  W.  Dietz,  Daniel  J.  Spohn , and  George  Merkel,  Status  of  the 
Tactical  Environment  Multiple  Systems  Evaluation  Program  (TEMSEP) , Harry 
Diamond  Laboratories  TM-77-23  (September  1977). 

‘*E.  R.  Parkinson,  SAPSC:  A Two-Dimensional  Close-in  EMP  Coupling 
Analysis  Code,  Science  Applications  Inc.  SAI-75-506-AQ  (October  1975). 
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JPLOT  If  equal  to  one,  the  transfer  function  amplitude  and 
phase  (if  used)  are  plotted;  otherwise,  not. 

IPHASE  If  equal  to  one,  phase  is  included  in  the  transfer 
function;  if  not,  amplitude  only. 

Card  No.  6;  FORMAT  (4(2E10.3));  (one  to  seven  cards) 

FTRANS(I)  Frequency  (Hz)  at  which  transfer  function  amplitude 
value  is  given. 

ATRANS(I)  Amplitude  value  (amps/volt)  of  measured  transfer 
function. 

Card  No.  7;  FORMAT  (IIP,  E10.3)  (used  only  if  IPHA5E=? 

NPHASE  Number  of  transfer  function  frequency  points  (25  or 
fewer) . 

TRETRD  Amount  of  time  (seconds)  by  which  the  input  pulse 
must  be  retarded  if  the  transfer  function  includes 
phase . 

Card  No.  8;  FORMAT  (4(2E10.3))  (one  to  seven  cards)  (used  only 
if  IPHASE=1) 

FPHASE(I)  Frequency  (Hz)  at  which  transfer  function  phase  value 
is  given. 

PHASE (I)  Phase  value  (degrees)  of  measured  transfer  function. 

4 . 2 Description  of  TIiANS2 

TRANS2  was  originally  written  to  operate  on  a CDC  6600  computer 
and  then  converted  to  run  on  an  IBM  370-series  computer.  Because  of  the 
short  word  length  of  the  IBM  computer,  it  was  felt  that  the  precision  of 
the  variables  and  functions  within  the  program  should  be  increased.  To 
avoid  the  necessity  of  ma)<ing  a substantial  number  of  changes  to  the 
program,  the  AUTODBL(DBL)  option,  which  automatically  doubles  the  pre- 
cision of  all  real  and  complex  variables  and  library  functions,  is  used. 
So  that  the  precision  of  the  plot  label  variables  is  not  doubled,  these 
variables  are  declared  to  be  integers. 

After  the  plot  labels  are  defined,  a number  of  arrays  are 
equivalenced  to  reduce  storage.  If  the  dimensions  of  any  of  these 
arrays  are  changed,  or  the  positions  of  the  arrays  within  the  program 
are  changed,  the  EQUIVALENCE  statements  must  be  altered  to  reflect 
these  changes. 


! 


I The  first  data  card,  containing  NDSETS  and  ISAME,  is  then  read. 

' NDSETS  is  the  number  of  data  sets  on  file  which  are  to  be  processed,  and 

! the  remainder  of  TRANS2  consists  of  a DO  loop  that  is  executed  NDSETS 

j times,  once  for  each  data  set. 

The  first  statement  inside  the  loop  calls  subroutine  READTP, 
which  reads  the  time  (T)  and  amplitude  (FT)  data  pairs  from  logical  unit 
No.  8.  If  the  number  of  data  pairs  (NPTS)  exceeds  300,  the  program  is 
terminated  because  the  dimensions  of  T and  FT  would  be  exceeded.  Next, 
the  input  parameters  are  read  and  printed,  unless  ISAME=1  and  other  than 
the  first  data  set  is  being  processed. 

NOFFR  data  points  are  then  removed  from  the  rear  of  the  pulse. 
This  feature  is  useful  if  the  oscilloscope  picture  is  accidentally 
digitized  beyond  the  region  of  interest.  The  pulse  is  inverted,  if 
desired,  ctnd  NOFF  extraneous  points  are  removed  from  the  front  of  the 
pulse. 

If  phase  is  included  in  the  transfer  function,  the  digitized 
pulse  is  now  retarded  by  TRETRD  seconds.  Since  it  takes  a finite  amount 
of  time  for  the  sensor  to  respond  to  the  input  signal  and  for  the  sensor 
output  to  pass  through  the  cables  and  the  balun  before  it  reaches  the 
oscilloscope,  the  time  frame  at  the  oscilloscope  is  retarded  from  the 
time  frame  at  the  sensor.  Because  use  of  the  phase  relationship  between 
the  sensor  input  and  the  cable  output  would  advance  the  output  forward 
into  the  sensor  time  frame,  the  digitized  pulse  must  be  retarded  so  that 
it  is  not  advanced  into  negative  time. 

The  amount  of  retardation  needed  for  a sensor-cable-balun 

system  can  be  measured  approximately  by  sending  a pulse  through  the 

system.  Unfortunately,  the  amount  of  retardation  depends  to  some  extent 

upon  the  frequency  content  of  the  pulse.  Because  of  this  dependence,  it 

may  be  necessary  to  make  a preliminary  run  of  TRANS2,  with  TRETRD 

slightly  larger  than  the  expected  value.  After  the  processed  pulse  is 

examined  to  determine  at  what  time  it  actually  starts,  TRETRD  can  be 

adjusted  accordingly,  and  TRANS2  can  be  run  again  with  the  correct  value 

of  TRETRD. 
s 

Next,  the  amplitudes  of  the  first  and  last  digitized  points  are 
set  equal  to  zero,  if  desired,  and  the  retarded  data  are  plotted  with 
the  DRAW4  plotting  routine,^  if  desired. 


^Thomas  V.  Noon, 
Libraries  ANAPAC  and 
1976) . 
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The  digitized  data  are  then  multiplied  by  sensor  response 
functions  which  are  not  included  in  the  transfer  function.  The 
following  list  describes  the  response  functions  presently  being  used. 

(1)  The  current  probe  output  is  multiplied  by  the  number  of 
resistors  (20)  in  the  banjo  to  provide  the  total  coupled  current. 

(2)  The  current  probe  output  is  simply  multiplied  by  1,  since 
the  total  current  is  measured  by  the  probe. 

(3)  The  current  probe  output  is  divided  by  the  effective  area 
of  the  coil  to  provide  the  Compton  current  density.  At  present,  an 
effective  coil  diameter  of  1.875  in.  (4.76  cm)  has  been  assumed  for  the 
current  probe  with  a center  opening  of  1.625  in.  (4.13  cm). 

(4)  Because  the  transfer  function  of  the  Moebius 
loop-cable-balun  system  cannot  be  measured  directly,  the  Moebius  loop  is 
replaced  by  a current  probe  to  ma)ce  the  transfer  function  measurement. 
If  this  transfer  function  is  multiplied  by  the  impedance  of  the  current 
probe,  the  result  is  the  transfer  function  of  the  cable-balun  system. 
Multiplying  this  result  by  the  sensor  response  of  the  Moebius  loop 
(1  weber/m^s)/(l. 14  x 10"^  V)  then  provides  the  value  of  B-dot.  At 
present,  the  impedance  of  the  current  probe  is  assumed  to  be  simply 
2.4  0,  although  it  is  actually  a function  of  frequency. 

(5)  The  current  probe  output  is  simply  multiplied  by  1,  since 
the  total  current  is  measured  by  the  probe . 

(6)  If  the  transfer  function  is  measured  with  a current  probe, 
the  output  of  the  one-turn  Rogowsici  coil  must  be  multiplied  by  the 
impedance  (2.4  S2)  of  the  current  probe.  The  output  of  the  one-turn 
Rogows)ci  coil  must  then  be  multiplied  by  the  sensor  response  of 
1.  X l0"^/volume  to  provide  the  Compton  current  density.  The  volume 
presently  being  used  is  5 x lo'^  m^ . 

Next,  the  frequency  spacing,  DF,  and  the  time  spacing,  DT,  are 
calculated.  The  frequency  spacing  is  the  interval  between  data  points 
after  the  data  are  transformed  into  the  frequency  dom.ain,  and  the  time 
spacing  is  the  interval  between  the  data  points  after  the  data  are 
transformed  bac]c  into  the  time  domain. 

The  data  are  then  linearly  interpolated  with  CLINTD,^  which 
also  extends  the  data  cirray  to  fulfill  the  requirement  that  the  maximum 
time  value,  TMAX,  equals  1/DF.  The  number  of  equispaced  time  points 

^Thomas  V.  Noon,  User's  Manual  for  the  Modular  Analysis-Package 
Libraries  ANAPAC  and  TRANL,  Harry  Diamond  Laboratories  TR-1782  (November 
1976) . 
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extending  across  the  range  of  the  original  data,  NTP,  is  then  found  so 
that  the  proper  time  interval  of  the  processed  data  can  be  plotted 
later.  The  equispaced  data  are  then  transformed  into  the  frequency 
domain  by  FFTD,  which  is  a double-precision  version  of  FFT. ^ 

The  transfer  function  amplitude  and  phase  information  is  avail- 
able at  only  a limited  number  of  frequencies.  To  provide  transfer  func- 
tion information  at  all  the  frequency  values  used  by  FFT,  a cubic  spline 
fit  to  the  data  is  performed.  Initially,  the  spline  coefficients  are 
calculated  by  the  subroutine  SPLICO,  which  uses  the  technique  discussed 
by  Pennington.^  The  coefficients  for  the  amplitude  and  the  phase  values 
of  the  transfer  function  are  stored  in  the  arrays  COEF  and  COEFP, 
respectively. 

Various  initializations  are  performed,  and  DO  loop  30  is  then 
executed  to  multiply  the  frequency-domain  representation  of  the  data  by 
the  transfer  function.  Since  half  aliasing^  is  used  by  the  Fourier 
transform,  only  the  first  half  of  the  data  array  plus  one  point  (NS21 
points)  is  multiplied  by  the  transfer  function. 

If  the  frequency  value  is  outside  the  range  of  the  spline  fit 
to  the  transfer  function  amplitude  or  phase  information,  the  subroutine 
ENDFIT  is  called.  During  the  first  execution  of  ENDFIT  for  einy  particu- 
lar fit  (KFIT=1) , fit  coefficients  are  calculated,  and  a calculation  of 
the  appropriate  transfer  function  value  is  made.  During  subsequent 
executions,  these  coefficients  are  not  calculated  again,  but  a calcula- 
tion of  the  transfer  function  is  made.  The  following  four  fits  are 
presently  being  used; 

(1)  The  transfer  function  amplitude  value  is  set  equal  to  1000 
if  the  frequency  is  below  that  of  the  spline  fit  information.  It  has 
been  assumed  that  the  user  will  supply  transfer  function  amplitude 
information  for  at  least  one  frequency  point  below  the  second  frequency 
value  (equal  to  the  frequency  spacing  1 x DF)  used  by  FFT,  so  that  the 
only  low-frequency  value  for  which  an  amplitude  value  is  needed  is  0 Hz. 
It  was  felt  that  some  large  value,  such  as  1000,  was  appropriate  for  the 
0-Hz  frequency  point,  since  the  balun  should  not  transmit  dc  signals. 

(2)  A decreasing  exponential  is  fit  to  the  high-frequency  end 
of  the  transfer  function  cunplitude  information.  Since  the  Tektronix  556 
oscilloscopes  used  at  AURORA  have  a high-frequency  3-dB  cutoff  at 
50  MHz,  it  was  felt  that  any  energy  present  in  the  digitized  data  at 

^Thomas  V.  Noon,  User's  Manual  for  the  Modular  Analysis-Package 
Libraries  ANAPAC  and  TRANL,  Harry  Diamond  Laboratories  TP-1782  (November 
1976)  . 

^Ralph  H.  Pennington , Introductory  Computer  and  Numerical 
Analysis  (2nd  ed.),  MacMillan,  Inc.  (1970). 
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frequencies  much  beyond  50  MHz  is  probably  due  to  the  process  of 
digitizing  the  data  and  should  be  filtered  out.  Therefore,  it  is 
suggested  that  the  last  value  of  the  measured  transfer  function  which 
the  user  provides  to  TRANS2  be  at  a frequency  of  50  MHz.  ENDFIT  will 
then  exponentially  decrease  the  value  of  the  transfer  function  cimplitude 
until  it  reaches  0 at  100  MHz.  Due  to  sampling  theorem 
considerations,®  the  value  of  FMAX  should  always  be  at  least  twice  the 
value  of  the  frequency  at  which  the  transfer  function  amplitude  drops  to 
0. 


(3)  The  transfer  function  phase  is  set  equal  to  180  deg  if  the 
frequency  is  below  that  of  the  spline  fit  information.  Again,  it  has 
been  assumed  that  the  user  will  supply  at  least  one  transfer  function 
value  below  a frequency  equal  to  1 x df.  Examination  of  the  output  of 
the  network  analyzer  indicates  that  a value  of  180  deg  at  0 Hz  is 
reasonable. 

(4)  The  transfer  function  phase  is  fit  to  a straight  line 
passing  through  the  last  two  points  if  the  frequency  is  above  that  of 
the  spline  fit  information.  Examination  of  the  output  of  the  network 
analyzer  indicates  that  this  fit  is  reasonable. 

If  the  frequency  value  is  within  the  range  of  the  transfer 
function  information  provided  by  the  user,  the  subroutine  SPLINE  is 
called  to  calculate  the  proper  value  of  the  transfer  function  using  the 
coefficients  calculated  by  SPLICO.  The  equations  used  are  presented  in 
Pennington's  book,®  but  the  search  technique  is  not  the  same  as  that 
presented  by  Pennington.  The  frequency  at  which  the  treinsfer  function 
is  needed  is  compared  with  increasing  values  of  the  frequency  intervals 
of  the  cubic  splines  until  the  proper  frequency  interval  is  found.  The 
transfer  function  is  calculated  using  the  spline  coefficients  for  the 
proper  interval,  and  a pointer  (KSTART)  is  used  to  "remember"  this 
interval.  During  the  following  call  of  SPLINE,  the  search  begins  with 
the  interval  designated  by  KSTART,  rather  than  with  the  first  interval. 
KSTART  is  updated  and,  during  each  subsequent  call  of  SPLINE,  the  first 
interval  used  during  the  search  is  the  interval  used  for  the  SPLINE  fit 
during  the  last  call.  This  process  works  because  the  values  of  the 
frequency  increase  with  each  call  of  SPLINE. 

After  the  correct  values  of  the  transfer  function  amplitude  cuid 
phase  values  are  calculated,  the  frequency-domain  representation  is 
multiplied  by  the  transfer  function,  given  by 


^Ralph  H.  Pennington,  Introductory  Computer  Methods  and  Numerical 
Analysis  (2nd  ed.) , MacMillan,  Inc.  (1970). 

^George  R.  Cooper,  and  Clare  D.  McGillem , Methods  of  Signal  and 
System  Analysis,  Holt,  Rinehart,  and  Winston,  Inc.  (1967). 
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TF  = Ae"^® 


where  A is  the  amplitude  (in  amps/volt)  and  0 is  the  phase  (in  radians) 
of  the  transfer  function. 

Before  the  end  of  DO  loop  30,  various  arrays  are  filled  for 
plotting.  FPLT  contains  the  frequency  points,  AMPPLT  contains  values  of 
the  transfer  function  amplitude,  PHIPLT  contains  values  of  the  transfer 
function  phase  (in  degrees) , and  FAMP  contains  the  frequency-domain 
representation  of  the  processed  data. 

If  the  user  wishes,  the  transfer  function  amplitude  and  phase 
values  and  the  frequency-domain  results  are  then  plotted  as  a function 
of  frequency.  If  ISAME=1,  the  transfer  function  is  plotted  only  once. 

Next,  FFTAD,  a double-precision  version  of  FFTA,^  is  called  to 
transform  the  frequency-domain  results  into  the  time  domain.  Because  an 
inverse  FFT  can  cause  an  overall  amplitude  shift,  the  amplitude  at  the 
first  time  point  (FTTl)  is  subtracted  from  all  the  other  values  of  the 
processed  data. 

The  processed  data  are  then  plotted.  If  ISENSR=1  and  IPLTC=1, 
coupled  current  predictions  made  by  the  computer  program  SAPSC**  are 
plotted  on  the  seime  graph  as  the  processed  data.  This  SAPSC  information 
is  read  from  logical  unit  No.  10.  Since  the  plot  array  (CURRNT)  is 
equivalenced  to  the  first  SAPSC  coupled-current  array  (CURNTl) , the 
current  predicted  for  measurement  position  1 is  used  unless  NCUR  is  not 
equal  to  1.  If  NCUR  is  not  equal  to  1,  the  plot  array  (CURRNT)  is  set 
equal  to  the  SAPSC  coupled-current  array  (CURNT2)  for  position  2. 

When  Moebius  loop  (ISENSR=4)  data  are  processed,  not  only  B-dot 
is  plotted,  but  also  the  tangential  H-field  is  plotted.  H is  calculated 
by  dividing  B by  the  permeability  of  free  space  (1.26  x 10"^  henry /m) , 
and  integrating  the  result  by  using  the  trapezoidal  rule.^ 

When  one-turn  Rogowski  (ISENSR=6)  coil  data  are  processed,  Isoth 
J and  J-dot  are  plotted.  J is  calculated  by  using  the  trapezoidal  rule 
to  integrate  the  J-dot  results. 


^Thomas  V.  Noon,  User's  Manual  for  the  Modular  Analysis-Package 
Libraries  ANAPAC  and  TRANL,  Harry  Diamond  Laboratories  TR-1782  (November 
1976) . 

**£.  R.  Parkinson,  SAPSC:  A Two-Dimensional  Close-in  BMP  Coupling 

Analysis  Code,  Science  Applications , Inc.  SAI-75-506-AQ  (October  1975). 

^Francis  Scheid,  Schaum's  Outline  of  Theory  and  Problems  of  Numerical 
Analysis,  McGraw-Hill  (1968). 
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Finally,  if  desired,  the  processed  data  are  written  to  logical 
unit  No.  9 so  that  they  can  be  cataloged  as  a disk  file  and  used  at  some 
future  time. 

4 . 3 TRANS2  Sample  Deck 

A sample  deck  needed  to  run  TRANS2  from  an  UPDTE  file  is  shown 
in  figure  4.  Information  concerning  the  lEBUPDTE  utility  has  been  given 
by  IBM.^  If  TRANS2  is  to  be  run  as  shown  in  figure  4,  at  least  one 
CHANGE  card  must  be  present  in  the  deck,  even  if  the  chcuige  is  merely  to 
insert  a comment  card  somewhere. 

The  data  set  name  (DSN)  associated  with  logical  unit  No.  8 
should  correspond  to  the  file  containing  the  data  to  be  processed.  The 
DSN  associated  with  logical  unit  No.  10  should  correspond  to  the  file 
created  by  the  computer  progrem  SAPSC.  If  SAPSC  output  is  not  to  be 
plotted  on  the  seune  graph  as  the  processed  data,  the  DD  cards  associated 
with  logical  unit  No.  10  should  be  removed  from  the  deck.  The  DSN 
associated  with  logical  unit  No.  9 should  be  the  user's  choice  of  the 

//JFD5CTR2  JOB  (HK3020 , , 5) , ' JFWD I ETZ ' ,CLASS=C ,MSGLEVEL=  (2  , ) ) 

//  EXEC  PGM=IEBUPDTE 

//SYSPRINT  DD  SYS0UT=A 

//SYSUTl  DD  DSNAME=HK302O.TRANS2,DISP=(OLD,PASS,CATLG) 

//SYSUT2  DD  DSNAME=6&NEWPL,DISP=(NEW.PASS) ,UN I T=V I 0 , DCB=* . SYSUTl , 

//  SPACE={TRK, (10,10)) 

//SYS  IN  DD  * 

./  CHANGE 

changes 
./  ENDUP 

//  EXEC  ANAF0RT,PARM.F0RT='AUT0DBL(DBL) ' ,F1=,0UT=X 

//FORT.SYSIN  DD  DSN=&6NEWPL , D I SP= (OLD , PASS) 

//G0.FT08F001  DD  DSN=HK3020 . CR35 , D I SP=SHR 

//GO.FT09F001  DD  DSN=HK3020 . TRANSOUT ,D I SP=(NEW,CATLG) ,UNIT=SYSDA, 

//  SPACE=(TRK, (1 ,1)) ,V0L=SER=USER02 
//GO.FTlOFOOl  DD  DSN=HK3020 . SAPSCOUT ,D I SP=SHR 

//GO.SYSIN  DD  * 

data  cards 


/* 

// 

Figure  4.  Sample  deck  for  TRANS2. 


^IBM  OS/VS  Utilities  Manual,  GC35-0005-4  (September  1976). 


name  of  a file  on  which  the  processed  data  are  to  be  stored.  If  the 
user  does  not  wish  to  retain  the  processed  data  on  a disk  file,  the  DD 
card  associated  with  logical  unit  No.  9 should  be  removed  from  the  deck. 

4.4  TRANS2  Sample  Plots 

Figures  5 through  9 contain  examples  of  some  of  the  plots  which 
can  be  produced  by  TRANS2. 


Figure  5.  Typical  plot  of  digitized  data  after  they  have  been  shifted. 


FDEOUENCY  (Hi) 


Figure  6.  Typical  plot  of  transfer  function  amplitude 


FREQUENCY  (Hz) 

Figure  7.  Typical  plot  of  transfer  function  phase. 
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5 . PROGRAM  PUNCH 


If  a permanent  record  of  the  processed  data  is  needed,  the  program 
PUNCH  can  be  executed.  A listing  of  PUNCH  is  presented  in  appendix  C, 
and  figure  10  shows  a sample  deck  needed  to  run  PUNCH. 


//JFDIAPUN  JOB  (HK3020) , ' JFWDIETZ' ,CLASS=A,MSGLEVEL=(2.0) 
//STEP)  EXEC  F0RTPCLG,PARM.F0RT=' NOSOURCE, NOMAP' 

//FORT. SYS  IN  DO  * 


source  deck 


/••• 

//GO.FTlOFOOl  DD  0ISP=SHR.0SN=HK3020.TRANS0UT 

// 

Figure  10.  Sample  deck  for  PUNCH. 

5 . 1 Description  of  PUNCH 

PUNCH  transforms  the  processed  data  file  created  by  TRANS2  into 
punched  cards.  The  first  punched  card  contains  the  title  punched  in 
20A4  format.  The  second  card  contains  the  number  of  data  pairs  punched 
in  110  format.  All  subsequent  cards  contain  the  time  and  amplitude  data 
pairs  punched  in  3{2E12.4)  format.  At  present,  'PUNCH  converts  the 
processed  data  from  double  to  single  precision  before  punching  them. 

5.2  PUNCH  Sample  Deck 

The  DSN  associated  with  logical  unit  No.  10  should  be  the  name 
of  the  file  containing  the  information  to  be  punched. 


6 . PROCEDURE  PURGE 

It  is  suggested  that,  to  save  storage  charges,  all  data  files  be 
deleted  after  the  data  have  been  processed.  This  deletion  can  be  accom- 
plished easily  with  the  PURGE*  procedure. 


*From  recommendations  by  Egon  Marx  and  Thomas  V.  Noon,  Harry  Diamond 
Labor a tories  (1976) . 
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7.  CONCLUSIONS 


The  program  TRANS2  has  been  used  successfully  to  tr^sform  digitized 
data  taken  at  AURORA  into  the  parameters  being  measured,  such  as  coupled 
currents,  Compton  currents,  etc.  When  phase  is  included  in  the  trcinsfer 
function,  the  user  may  initially  have  some  difficulty  choosing  thfe 
proper  time  shift  to  use,  but  as  the  user  becomes  more  familiar  with  the 
program  and  the  data  being  processed,  this  choice  should  become  less 
difficult.  I would  suggest ^that  any  potential  user  of  TRANS2  have  some 
knowledge  of  Fourier  transforms  ^ and  the  sampling  theorem  before 
attempting  to  use  this  program.  A poor  choice  of  some  input  patrameters, 
such  as  NSTAR  and  FMAX,  could  result  in  the  data  being  improperly 
processed. 


LITERATURE  CITED 


(1)  Thomas  V.  Noon,  User's  Manual  for  the  Modular  Analysis-Pack&gc 
Libraries  ANAPAC  and  TRANL,  Harry  Diamond  Laboratories  TR-1782 
(November  1976)  . *'  . 

J2)  IBM  OS/VS  Utilities  Manual,  GC35-0005-4  (September  1976). 

(3)  John  F.  W.  Dietz,  Daniel  j.  Spohn,  and  George  Merkel,  Status  of  the 
Tactical  Environment  Multiple  Systems  Evaluation  Program  (TEMSEP) , 

• Harry  Dicimond  Laboratories  TM-77-23  (September  1977). 

(4)  E.  R.  Parkinson,  SAPSC s A Two-Dimensional  Close-in  EMP  Coupling 
Analysis  Code,  Science  Applications  ■ Inc.  SAI-75-506-AQ 
(October  1975) . 

(5)  Ralpl.  H.  Pennington,  Introductory  Computer  Methods  and  Numerical 
Analysis  (2nd  ed.),  MacMillan,  Inc.  (1970). 

(6)  George  R.  Cooper,  and  Clare  D.  McGillem,  Methods  of  Signal  and 
System  Analysis,  Holt,*  Rinehart,  and  Winston,  Inc.  (1967). 

(7)  Francis  Scheid,  Schaum's  Outline  of  Theory  and  Problems  of  Numerical 
Analysis,  McGraw-Hill  (1968). 


22 


APPENDIX  A.--A  LISTING  OF  THE  PROGRAM  OFILE 
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